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Abstract 

We calculate all the four-photon helicity amplitudes at the one-loop level in a massive 
theory using multiple-cut methods. The amplitudes are derived in scalar QED, QED and 
QED^=^ theories. We will see the origin of rational terms. We extend the calculation to 
the simplest six-photon helicity amplitude where all photons have the same helicity. 
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1 Introduction 



The light-by-light scattering is a good laboratory to find efficient methods to compute a mas- 
sive loop multi-leg amplitude, because gauge invariance and IR/UV finiteness lead to enormous 
cancellations when we sum all the Feynman diagram. The first calculation of the four-photon 
amplitude in massive QED was done by Karplus and al. [Tj. They straightforwardly calcu- 
lated each Feynman diagram. Then B. de Tollis [2j computed the four-photon amplitude with 
Cutkosky rules. Ten years ago, Bern and Morgan calculated the four-gluon helicity amplitudes 
with a massive loop [3]. They used the two-cut unitarity methods in n dimensions. Recently, 
Bern et al. calculated the two loop corrections QCD and QED to light-by-light scattering by 
fermion loops in the ultrarelativistic limit [4j . Then Binoth and al. [5] calculated the four-photon 
helicity amplitudes in QED, scalar QED, supersymmetric QED-'^^^ and QED-'^^^ in massless 
theory. They used a n-dimensional projection method. Two years ago, Brandhuber and al. 
calculated the four-gluon one loop helicity amplitudes with generalized unitarity cuts [6]. 

Here we want to compute the four-photon amplitude, in massive QED, with a recent method: 
the generalized unitarity cuts. Even at the energy of LHC, 14 TeV, the mass of heavy quark t can 
have a significant effect. This article aims at explaining how to use new unitarity-cut methods in 
the case of massive theories. Here we apply this new technology for a 2 ^ 2 process. However, 
the future project is to calculate easily some 2 — > 4 QCD processes, like gg — > tttt with heavy 
quarks, present in the background of LHC. The knowledge of the background is very important 
to detect new particles like Higgs. 

In this paper, we calculate the four-photon amplitude at one loop order in three massive QED 
theories: scalar QED, QED and supersymmetric QED-'^^^. Our result agree with [4il5]. We call 
j^^caiar (^respectively j^^p^"-"^ and A^^^ ) the four-photon amplitude in scalar QED (respectively 
QED and supersymmetric QED^'^^^). We obtain very compact expressions contrary to Karplus 
and we can deduce easily the origin of the rational terms. Actually, we can link easily these three 
amplitudes with a supersymmetric decomposition. All diagrams of the four-photon amplitude 
in QED have the same pattern: four external photons coupled to a fermion loop. However, 
using the fact that degrees of freedom for internal lines can be added and subtracted [Hi [10], we 
write the internal fermion as a supersymmetric contribution and a scalar: 

/ = _2s + (/ + 2s) Ai""'" = -2 Af"^"^ + A-^=\ (1) 

This formula is true for massless and massive theories. Moreover we point out that this formula 
([1]) is true for any number of photons. 

We calculate all Af^"'^'^^ helicity amplitudes with generalized unitarity cuts in sections l3.1l 3.2| 3.31 
Then, we use extensively the supersymmetric decomposition ([T]) to calculate the QED ampli- 
tude A^^^^^°^ in section HI and the amplitude A^^^ in section [6| Indeed, if we write the 
four-photon amplitude in QED, the supersymmetric decomposition imposes the pattern of the 
supersymmetric amplitude A^^'^ function of magnetic moments. We don't need to calculate 
all the supersymmetric diagrams. After, we discuss about the origin of the rational terms and 
the different cut-techniques in section [5j Finally we derive the most simple helicity amplitude 
of the six-photon amplitudes in Section [71 

But first, some notations and explanations on generalized unitary cuts will be introduced. 
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2 Notations and explanations 



2.1 The structure of the amphtude 71 + 72 + 73 + 74 ^ 



We study the process 71 + 72 + 73 + 74 ^ 0. The momenta of the ingoing photon called i 
is pf. In this paper we suppose that all the photons are on-shell, so we have the first relation 
Vi G [1..4], p2 = 0. Dia grams at tree order are impossible and the first non vanishing order is 
one-loop order. 

The QED Lagrangian contains one vertex, whereas the scalar QED Lagrangian gives us two 
vertices. We recall them in Appendix [Aj Therefore in QED we have six one-loop diagrams 
whereas in scalar QED we have twenty one one-loop diagrams. As we have only four external 
particles entering in the loop, the power counting tells us that individual diagrams are UV 
divergent. So we regularized the divergence in calculating loops in n = 4 — 2e dimensions. 
Nevertheless, thanks to gauge invariance, the sum over all diagrams has no UV divergences. So 
we have to observe compensations. Moreover, as we suppose that the four photons are on-shell, 
therefore we decide to decompose the amplitudes j^^'^'^"^ ^ j^scaiar A^^^ on a basis of master 
integrals in n and n + 2 dimensions: 

= ^ (04/4 + ftj/g (Imass) • + Ci/2 (Imass) •) + rational terms, (2) 

iGl,2,3,4 

where I^^"^ is the no external mass box in n -|- 2 dimensions, I3 (Imass) is the three-point 
function with one external mass and I2 (Imass) is the two-point function with one external 
mass. The exact definition of I^'^'^, /g (Imass) and I2 can be found, for example in [3l EJ El [8] ; 
nevertheless to be self consistent we recall them in Appendix O The interest of this basis 
of master integrals, is to separate IR/UV, rational and analytic terms. But this basis is not 
unique. I^'^'^ has an analytic structure with polylogarithms whereas the IR divergences, in 
massless theories are carried by the function I^{lm) and the UV one, in massless and massive 
theories by the function . Each diagram is UV divergent so each diagram has scalar bubbles 
I2, but the amplitude (sum of diagrams) is UV finite so we expect to have compensations 
between the different scalar bubbles to eliminate the divergences. 

The amplitude is totally defined by the coefficients aj,6j,Cj and the rational terms. To 
calculate all these coefficients, we use the spinor formalism and the method of the helicity 
amplitudes. 



2.2 Spinor formahsm and hehcity amphtudes 

We use the spinor helicity formalism developed in For the spinorial product, we 

introduce the following notation: 

{Pa - \Pb+) = {ah) (3) 
(pa + \Ph-) = [ah] (4) 

iPa - \MPc-) = {ahc) = [cba] = {pc + |^b|Pa+) = [PcPb]{PbPa) (5) 
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{Pa + \ AT^c\Pd-) = [abed] = -[dcha] = -{pa + \i>ci)h\Va-)- (6) 

Moreover we introduce the classical Mandelstam variables : 

s = S12 = (12) [21] t = S14 = (14) [41] u = sxz = (13) [31]. (7) 

All coefficients a^, 6i, Cj and rational terms are described as products of spinor and Mandelstam 
variables. 

We calculate all helicity amplitudes. A photon has two helicity states a = it. The amplitude 
is the sum of all helicity states. 

^4 = ^ ^4 (o-i, 0-2, 0-3,0-4) . (8) 

(Ti=±,o-2=lt,(T3=±,(T4=± 

As we have two helicity states per external photons and four external photons, so the amplitude 
is the sum of 2^ = 16 helicity states. However they are not all independent. In fact we have 

only three independent helicity states ^4(+ + ++),^4( — V ++) and ^4( Others are 

obtained by permutations and parity. 

It is necessary to introduce the polarisation vectors of the external photons. The advantage of 
the helicity amplitudes is that we can express easily all the photon-polarisation vectors with 
spinors. Those expressions come from [11]: 

^ ^(m) ^ \/2[lr]' ^ ' 

where ji?) and |r) are two arbitrary light-like vectors. Before giving results on A^'^"'^"'^, we give 
explanations and notations on regularization, propagators and integrals. 



2.3 The regularization scheme and notation of integrals 



All diagrams have the same pattern and all are UV divergent. To regularize, we put the 
loop momentum Q'^ in n = 4 — 2e dimensions. Capital letters describe vectors in n dimensions 
whereas small letters describe vectors in four dimensions. So we decompose the loop momentum 

= qf^ fj^f^ with q^^ the four-dimensional part and fi^ the — 2e part. The four-dimensional 
space and the — 2e-dimensional space are orthogonal so one has: = — fi^. And we operate 
in the "four-dimensional helicity scheme" , in which all external momenta are in four dimensions. 
In n dimensions, we take the prescription where the propagators are defined by: 

fermion i— 75 ^ — = t-^ ^ ^ — = i — —5 — , (10) 

QZ _ _|_ gZ _ _ j^Z _|_ jJZ 

.1 _ . 1 _ . 1 

scalar ^-Tv) o ^ — o 5 ^ — ^ . (H) 

QZ _ j^Z l\ gZ _ ^Z _ j^Z _|_ JJZ 

The formulas (jlOl lip show that the running particle seems to have a mass: + We 
will observe this phenomena in the analytic expression of the amplitude A4. The generalized 
unitarity-cut imposes the calculation of some trees in n dimensions, where the running particle 
enters with a mass m? + /i^. The n-dimensional calculation introduces some integrals with 
powers of in the numerator; those integrals are called extra-dimension integrals. We express. 
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in Appendix ID. 21 those integrals in terms of higher dimension loop scalar integrals. But here 
we give the definitions of some scalar integrals and extra-dimension scalar integrals, used in this 
paper [3]: 

= T^h"^^' 

We explain in Section [5.21 the origin of those extra scalar integrals. In our convention, we ignore 
a factor K = z(47r)~"/'^ in front of each scalar integral of the amplitude, and we reintroduce it, 
in the final result. In the following, we often add arguments to the integral to denote the order 
of photons entering the loop. We give the definition of those arguments: 

Pc. 

12 (abed) = l2{Sab,Sad) = 




I3M 



I^M = T-{sa,) = I >Q<1 (15) 

In Appendix [Cl we give some analytic expressions in terms of polylogarithms of those massive 
and massless scalar integrals. 

We decide to express the amplitude ^/e'"™«o?i/scaZar/A^=i combination of master integrals 
([2]). To calculate the different coefficients in front of each master integral, we use the generalized 
unitarity-cut in n dimensions, which we recall in the next subsection. 



2.4 Generalized unitarity-cuts 



The unitarity cut method come from the Cutkosky rules [T2|. In the last ten years, there 
has been an intense development around the unitarity-cuts and several generalizations have 
been made. The first generalization of those rules is that we can cut not only two propagators 
but also three [131 E! 01 four propagators [151 HSl E] • But we will see that the more we cut 
propagators, the more we loose information. However, the more we cut propagators, the more 
the amplitude is simple to calculate. The second generalization is to evaluate the loop integral 
in n = 4 — 2e dimensions [3], which was improved in [TF]. The four-dimension cuts are very 
efficient to calculate coefficients in front of structures but the extra-dimensional cuts are powerful 
to calculate the rational terms. We will see a link between extra-dimensions and rational terms 
[191 120j . But the generalization of this link is not so obvious for a general amplitude. Another 
extension, which has recently be done, is the generalization of the unitarity cuts to massive 
theories. In [3], we find the calculation of four-gluon amplitudes at one loop in massive theory 
with the two-cut techniques. And recently, the unitarity-cut techniques in massive theories was 
generalized to three and four-cut techniques [2Tj. Finally, a few months ago, Papadopoulos 
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and al. gave a general method to calculate each coefficient in front of the master integrals [22| . 
extended by Forde [23l[2i]. 



The Cutkosky rules [12] require to consider an invariant or a channel, constituted of several 
consecutive legs. Consider a loop amplitude, called A. One first computes the discontinuities 
across branch cuts (imaginary parts) by evaluating a phase space integral. The imaginary part 
of the amplitude is the sum over all the discontinuities: 

2 Im(^r'"' (+ + ++)) = Yl Disc^A (16) 

iGchannel 

The discontinuity in a channel can be computed by replacing the propagators separating the 
set of legs by delta functions: 

^ ^ 2. 6M (Df) . (17) 

i 

The real part is then reconstructed via a dispersion relation. The existence of a linear com- 
bination of scalar integrals allows us to avoid this reconstruction explicitly. We perform the 
cut calculation a bit differently. Consider the loop amplitude, with a cut in the channel: si2 
(FigH]). We compute the discontinuity Disc^jj {A). It is convenient to replace the phase-space 
integral with an unrestricted loop momentum integral which has the correct branch cuts p9|) . 
In this integral ()19p . the tree-amplitudes are kept on-shell. Then we decompose the discon- 
tinuity DisCs^2 (^) ^ * linear combination of scalar cut-integrals in this channel (j20p . The 
reconstruction of the amplitude is hidden in the rebuilding of the scalar integrals. 



pi 
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|P4 
P3 



Fi gur6 1: Fermion loop with a cut in the channel si2. 



Disc,,, (A) = {2TTfj ^ 6^+\Dl) Atree{2) 5^+\DI) Atreeil) 



I 



Atreei'2) 7^^tree(l) 



(27r)"Di^ 



S12 



E^^^"i^i2 = E^^Disc,,^ (If). 



(18) 
(19) 
(20) 



From this, we can write: 



A = ^air + A. 



(21) 



In the term A, there is a combination of scalar integrals which cannot have a cut in the chan- 
nel Si2. To obtain the coefficient in front of all scalar integrals, we have to consider cutting 
amplitudes in all channels (si2, sia, S14). In the following, we note: 



2^5(+) 



5. 



(22) 



We introduce a notation to label the number of cuts and the channels. "Disc2,s" means cutting 
the two internal lines in channel "s". Cutting a third internal in all possible ways leads to 
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"Discs, s", while cutting the two remaining hnes leads to "Disc4", in which there is no need to 
specify the channel. We denote DIscat = X^iechannei -^i^'^A^.Sj, = 2,3. 

We are going to calculate all the helicity amplitudes with two, three and four-cut techniques 
in n = 4 — 2e dimensions for massive theories, and then we compare all those techniques. Now, 
as we have given the definition of all objects used in this paper, we are going to calculate the 
first helicity amplitude 2"*", 3"*", 4+) in the next section. 



3 ji^scaLar ^elicity amplitudes 



3.1 Af'^'"^(l+,2+,3+,4+) helicity amplitude 
3.1.1 Four-cut technique 

The four-cut technique [T5l[16l[T7] says that we cut all four propagators D^, i = [1..4]. We 
have: 



Disc4(Ar"^'-(++++)) = ^ 



a(l,2,3,4) ^ ^p. 

We define the loop momenta of propagators as Qi = Qi^i + pi and Qo = Qi- Using Feynman 
rules, we compute the discontinuity Disc (^|'^"'"^(-|- -|- -|-+)): 



Disc4 (^r'"'(+ + ++)) = j 



(i?3) {RA) 



-5 {Dt) 6 {Di) 5 {Di) 6 {Di) . (23) 



As the four-dimensional and the — 2e-dimensional spaces are orthogonal, therefore, the spinor 
product can be simplified: {RQ22) = {Rq2'^) + (i?/^2) = {Rq22). In the following, we do this 
simplification each time it is possible and then we use the fact that all propagators are on-shell 
to simplify eq. (f23|) . We use the first on-shell tree computed in Appendix [Fl This tree has 
only two photons, so we split the integrand into two groups of photons {pi,P2) and (^3,^4), 
and we apply the formula ()F.45|) for the two groups of photons. The discontinuity ()23p directly 
becomes: 



Disc4 {Ar'^^i+ + ++)) = (eV2)^ Y / d^'Q + ^'f ^ (^?) ^ m 5 {dI) 6 

(24) 



o-(2,3,4) 

(^^)' Y 7liSDisc4(K4"(1234)). (26) 

o-(2,3,4) ^ ' 
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3.1.2 Three-cut technique 



The three-cut technique imposes three propagators on-shell. So in the case of the four-photon 
amplitudes, we have four branch cuts per four-point diagram, only one branch cut per three- 
point diagram and zero branch cut per two-point diagram. But actually, as the four photons are 
on-shell, we have only two invariants per four-point diagrams, and therefore only two indepen- 
dent branch cuts. So we divide the result by two. Moreover we collect diagrams to construct set 
of gauge invariant trees, so we divide again per two. The discontinuity Discs (^f + ++)) 
is: 

Discs (^r'"'(++++; 




In the first group, for example, Df is the propagator between the two photons pi and p2 but 
we should not forget the diagram with the four-point vertex. There are three diagrams in each 
group. Thanks to permutations, all groups of cut-diagrams are the same, so the discontinuity 
becomes: 

Discs {Ar'^^{+ + ++)) =\ J}f"' • (27) 

^(1,2,3,4) \^ 

So using the Feynman rules, the discontinuity Discs (^4'^"''^'"(+ + ++)) is: 

<t(1,2,3,4) \(t{1,2) ^ ' 1 \ / 

(28) 

We use the expression of the on-shell tree ()F.45h for the second group (^3,^4). For the first 
group of photons (^1,^2)) the propagators around this group are on shell, but the propagator 
joining the two photons are not on-shell, so we use (|F.49P : 

<t{1,2) \ / 1 \ / \ / ^(1^2) 1 

Inserting (fR45]l and (l29|) in (l28|), the discontinuity Discs (Af + ++)) becomes: 

(t(1, 2,3,4) ^ ' •' o-(l,2) 



Discs (^r'-(+ + ++)) = ^ ^ E Uli) / i^' + E i^l) ' (^i) ^ (^l) 

(30) 



(^^)' E tJI^ (Discs,... (1234))+ Discs,.,, (i^r(1234))) 

<t{2,3,4) ^ ' 

(31) 
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= {eV2f 7iS°i^^3(i^r(1234)). (32) 

<t(2,3,4) ^ ' 

In the last step we have gathered the two branch cuts of the scalar integrals K"^ to rebuild the 
entire discontinuity. Indeed the computation could be more simple. Consider a discontinuity 
with several branch cuts. The self consistency of the unitarity says that we find the same 
coefficients in front of all branch cuts of each scalar integral of this discontinuity. Actually, we 
need to calculate only the coefficient in front of one branch cut of each scalar integral of the 
discontinuity. 



3.1.3 Two-cut technique 



This time, only two propagators are on-shell. As we have on-shell photons, we have only 
two channels per diagram. Therefore we have two branch cuts and the imaginary part of 
the amplitude is the sum of the discontinuity of the two branch cuts. So, for each diagram, 
we decide to cut the propagators Z)|,-D| and then we cut the propagators D\,D^. Finally 
regrouping diagrams with the same cuts to make trees, we obtain: 



2 Im (Ar'"''(+ + ++)) = Disc2 (^r''^'^(+ + ++)) = \y1 



Pi Dj pj p, Dj m 

D! |,l D| + D! -^-Dl 

'^(2,3,4) \p, £,2 t,p^ p,, ^ ^2 4>, 



(33) 



The numerical factor | comes from the fact that we have gathered diagrams to create gauge 
invariant trees. We point out, thanks to permutations, the two groups of cut-integrals are the 
same. Therefore the discontinuity Disc2 (^|™'"^(-|- + ++)) is written: 

<t{2,3,4) \f7{l,2) ^ ' 1 \ ' 



* 



(34) 

Using (lF.49p . to simplify the two trees in ([Ml), the discontinuity Disc2 (^f'"'"(+ + ++)) be- 
comes: 

Disc2 (^r'-(+ + ++)) = E iiS / i^' + -^)' E E (^2) ^ (^4) 

(t(2,3,4) ^ ' '' (t(1,2) 1 cr(3,4) 3 



<7{2,3,4) 



(35) 

2(e^/2)4 E |||^Disc2,.,,(i^r(1234)). (36) 



We have one branch cut of the scalar integrals. But we want to reconstruct the entire dis- 
continuity of the scalar integral K2- Using the conservation of energy-momentum of external 
moments, is invariant by permutation. So we can split the amplitude in two equal terms. 
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We transform, thanks to permutations, one of this term to obtain the second branch cut, and 
the discontinuity of the scalar integral appears: 

Disc2 (^r'"X+ + ++)) = (eV2)4 Yl TllS i^^^^2,s,, {K2 (1234)) + Disc^,.,, {K2 (1234))) 

<7(2,3,4) ^ ' 

(37) 

= (^^)' E 7i!S°i^^2(Kr(1234)). (38) 

<t{2,3,4) ^ ' 

With this last discontinuity we reconstruct directly the helicity amplitude A|'^'^'"''(+ + ++) 
by transform the discontinuity of each scalar integral to the scalar integral. We multiply the 
result by the factor K = z(47r)~"/^: 

^ ' o-(2,3,4) ^ ' o-(2,3,4) ^ ' 

where a = e^/47r. The reconstruction of scalar integrals is automatic, we don't need dispersive 
relation. There is full agreement with [3]. 



3.2 A|™'''''(l-,2+,3+,4+) helicity amplitude 
3.2.1 Four-cut technique 

The four-cut technique sets the four propagators D?, i = [1..4] on-shell: 

Disc4 [Ar'^^^i- + ++)) = \ ^:if ■ 

a(l,2,3,4) „i \t 

So the discontinuity Disc4 (Af '"''(- + ++)) is: 

Disc. (Ar'-{- + ++)) = £ / <||> (|ll> ^ («;) 5 8 («|) i , 

(40) 

We split the integrand in two groups of photons: (pi,P2) and (^3,^4). As the helicity of the two 
photons of the first group are different, therefore we use (IF.56P in the limit of all propagators 
on-shell: 

T&rilr = .j:Sf-o» = -p^((1..232) + („^ + ,„^) 12311), (41) 
and for the second group of photons, we use the relation (|F.45p . The discontinuity (j40p becomes: 
Disc4(Ar''^^(- + ++)) = {eV2f Y ^^ld'Vif^' + m'fs{Dl)6iDl)6iDl)6{Dl) 

+ («^)' E T^^Mi) I ^"^(1^2232) + m') 5 {Dj) 5 [dI) 5 {dI) 5 [dI) . 

(7(2,3,4) ^ ' 

(42) 
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Here, to complete the computation, we have two ways. The first way is to use the integration 
formula of the tensor integral (|D.33p . The second way is to use the fact that we have four 
on-shell propagators, which gives us four conditions. Those conditions are sufficient to define 
exactly the loop momentum. We explain this calculation in Appendix El We find: 

Disc4(^r'''^(- + ++)) = ieV2)^ E ^j^^'^'^'i^^ + il^^)- (43) 



3.2.2 Three-cut technique 



We have a photon with a negative helicity. So not to break the helicity symmetry, we gather 
cut-diagrams in three groups rather than two or four, and we multiply by a factor 1/2. The 
conservation of the symmetry allow us to rebuild easily the discontinuities with the branch cuts. 
In terms of cut-diagrams, the discontinuity Discs (Af + ++)) is written: 

Discs (^r'"'(- + ++)^ 



r'vV r'vV''' I'^^v 

<7(2,3,4) \p} rk >J Pt Dl Pi Dl ^Pi 



So, using the Feynman rules, the discontinuity Discs (^4'^"'"^( — h ++)) is: 

Discs (^r^'^'^ (-+++)) = \ 



(7(2,3,4) 

y(T(3,4) y 

We calculate the first tree - \i}-\{^2) + Scr(i,2) "^My^ ' imposing \r) = |2), \R) = |1), and 

using the formula ()F.56p in the limit of the propagator Dg on-shell, we obtain: 

.tb) t^"] ^1 ^^^^ — 




(1^2232) + (^2 + ^2) [231] (1^2312) + (/^^ + m^) [231] 



+ 



(231) ^ Dl D[' J ■ 

(45) 

To evaluate the on-shell trees with the same helicities, we use ()F.45| F.49p . For the two last 
trees ^"j^^^^ ^^[TTJ^ ^TTrT ^^2^^ 1 ^^^i use again the formula (|F.56p , but with very few spinor 
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manipulations we obtain: 



{Rq^^) [rgil] ^ t^'^^l] .2^ 2^[421] .... 

(i?4) [Ir] = (-^^^ + (42Ty - + ) (421) ' ^^'^ 

[Ir] (i?2) = + (23Ty - + ) ^'^^ 

We gather (|45l 461 47p in the discontinuity (j44|) . We obtain some linear tensor integrals, which 
we integrate with (ID.33| D.34j ). Some three-point linear-tensor integrals appear, but many of 
them are zero. Consider the three-point function with one external mass S23: H{s2'i){(li)- 
This triangle can be expressed as a linear combination of its two on-shell legs (Appendix IC]) : 
-^3('S23)('7f ) = A p'^ + B P4. As the momenta pi and p4 are light-like vectors, therefore the 
linear tensor integral /3(s23)((lq'f 4)) is zero. After integrating, we have a linear combination of 
discontinuities of four-point, three-point and two-point scalar cut-integrals. But some of those 
discontinuities are spurious. Consider the invariant Sij. Sij is a channel of the one external mass 
triangle 13(3 ki) if the external mass Ski is equal to the channel Ski = Sij. And there is the same 
argument for the bubbles. Integrals which don't respect this condition are spurious and we drop 
them. They appeared because we lifted a cut-integral to a Feynman integral. For the rational 
terms there is no problem. The reduction of the Feynman's integrals gives some extra-scalar 
integrals. The development of those scalar integrals in function of e create the rational terms. 
So we keep only the extra-scalar integrals verifying cuts. Keeping only those integrals with cuts, 
we rebuild the discontinuity. The discontinuity of four point scalar integrals need two branch 
cuts to be rebuilt, whereas, only one branch cut is sufficient to rebuild a three or two point 
scalar integrals. So, we obtain: 



Disc3 [Ar'"-{- + ++)) = E H^Di-'^a (k4"(1234) + + ^JSU) - 



<7(2,3,4) 



(48) 



o-(2,3,4) 



Pointing out that ^^^j ^^^^i) = {23) {421) ' ^^^^ gather discontinuities. We symmetrize the coeffi- 
cient in front of the three-point extra-dimension integral J3 . The discontinuity Discs [A^'^"'^°'^[ — h -|-+)) 
becomes: 

Discs (^r'^'^(- + ++)) = {eV2)' Yl S§iyDiscs fe(1234) + £ Jr(1234) 

cr(2,3,4) \ '\ ' V 



3.2.3 Two-cut technique 
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Each diagram has to be cut in the two channels, corresponding to the two branch cuts. In 
terms of cut-diagrams, the discontinuity Disc2 (^|'^'^'°''( — h ++)) is: 



Disc2 (^r^'^'^ (- + ++)) = i Yl i ""^li^}' + '''^S^''' 



<7{2,3,4) \pt ^4 \t Pi ^^\t 

If we do permutations, we can see easily that all cut-diagrams are doubled. We gather them 
and the discontinuity Disc2 [Af^'^'-°-^{ — h ++)) is: 

Disc. (^r'"'(- + = 5 E (-^»)'/'''«f-ra7iT+ ^ ' 

<7(2,3,4) \ 



[lr](i22) J^^^ [Ir] Df {R2) 



* 



(50) 

We use (|F.56p and (|F.49p to calculate the trees of the discontinuity. We obtain some tensor 
triangles, which we integrate with the formulas (|D.33| D.34p . We keep only cut-integrals which 
are not spurious, and we rebuild the discontinuities. Finally, the discontinuity ([50]) is: 

Disc2(^r'"'^(- + ++))= {eV2)' S§iyDi^^2 fc(1234) + 1^ Jr(1234) 

o-{2,3,4) \ '\ ' V 

-(^)^?<'))' 

The two-cut technique gives all information to reconstruct the helicity amplitude. We find 
straightforwardly: 

„2 I +2 I „,2^ 



(52) 



o-(2,3,4) 



3.3 Af^'^'Xl ,2 ,3+,4+) helicity amplitude 

This helicity amplitude is usually called the MHV (Maximal Helicity Violating) amplitude. 

3.3.1 Four-cut technique 

One of the difficulty of this helicity amplitude, is that we have two kinds of topologies of 
helicities. The helicities are either alternate or they are paired as shown eq. (I53p . We group 
diagrams according the topology of helicities and the four propagators Df, i = [1..4] are on-shell 
so the MHV discontinuity is: 

Pi Dl p+ PF -03 Pt 

Disc4(Ar»'-(- -++))= y: E "'^:f^ + E • (53) 

a(l,2) <7(3,4) p- %j <t{1,2) p. 
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We obtain: 



Disc4(^r'"' (--++; 



(t(1,2) (t{3,4) 



a(l,2)- 

= /1 + /2. (55) 

We have split the two topologies into two integrals: Ii and l2- Applying two times the tree 
formulas (1F.45| F.46p . the first topology Ii is directly: 

h = {eV2f E Tll|S^i^^4(^4(1234)). (56) 

a(l,2) <t{3,4) ^ J \ ' 

For the second topology /2, instead of using (lF.56p . we gather photon with the same helicity 
and we reduce directly: 

[rq^l] [rq^2] _ {lqiq^2) _ (12) (/i^ + m^) + (14^32) ^^^^ 



[Ir] [2r] [12] [12] 

{Rq2^) {RqA^) ^ [3g2g44] ^ [34] (/i^ + m^) + [3^224] 
(i?3) {RA) ~ (34) ~ (34) 



(58) 



Therefore with (|57| 58|) . I2 becomes a sum of four terms which we develop and integrate with 
the formulas ()D.33| D.34|l . Here as we have four cuts, all triangles and bubbles are spurious, 
because they don't have four cuts. So I2 is spelt: 

4(12)[34] [^._ , 2ut^._ , uH^ 



h = 2(g^2) [\2] (34) V^'''^' (^4"(1324)) + — Disc4 (7^(1324)) + ^Disc4 (/4"(1324)) 

(59) 

The amplitude contains a four-point scalar integral in n dimensions I4 . This integral (I4) in 
a massless theory has IR divergences. Each diagram of the four-photon amplitudes has no IR 
divergence. So those divergences should be compensated by other divergent integrals like three- 
point scalar integrals. If we have three-point integrals in massless theory, we have probably the 
same in a massive theory. To simplify the problem it is better to transform the n-dimensional 
four-point scalar integral into a, {n + 2)-dimensional integral, which is no longer IR divergent. 
This transformation {l^ I^'^'^) is given by the formula (ID.40p . Keeping only integrals with 
four cuts, we have: 



I2 = 2(e^/2)4|^|M (^Disc4 (i^r(1324)) - ^Disc4 (1^+2(1324))) . (60) 

So as the discontinuity Disc4 (^^^a/ar^ 1"+)) is the addition of the integral Ji, given in (I56p 

and the integral I2, given in (I60p . therefore we obtain: 

Disc4 (^-'^'"X- - ++)) = (e^/2)^j||M|-^Disc4(/r+^(1324))+ J] Disc4 (Kr(1234)) 



ct(2,3,4) 

(61) 
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3.3.2 Three-cut technique 



The discontinuity, after grouping together diagrams with the same cuts is: 

1 



Discs A 



\ scalar , 



+ 



+)) =iE E 

cr(l,2) o-(3,4) 



\ P2 Di 



E E E '■') 

o-(l,2) (t(1,3) o-(2,4) * ^ ' 



"3 Pi Dl pY 



P2 P| "5 P2 / 



Using Feynman rules, the discontinuity is: 

{-iV2e)^ 



Discs (^^™'''^(- -++; 



(t(1,2) o-{3,4) 



<7(1,2) 



r(l,2) cr(3,4) <t(3,4) ^ \ ' 1 \ / / L J L J 

[rS]{Rl) [rqil] i {Rq23) \ [rqs2] {Rq^A) 

[lr](i?3) ^^^^ [It] Dj {R3) j [2r] {R4) 



a{l,2) <7(3,4) 

+ E E E /''"<3 

(7(1,2) (t{1,3) cr(2,4) 



5 (i^i) 6 {Dl) 5 {Dj) 



C7(l,2) <7(1,3) (7(2,4) 



(7(1,4) 



(62) 



We are not going to develop all the computation because there is no difficulty and all trees 
have already been calculated in this paper. It remains some tensor integrals, which are reduced 
with the formulas (ID.33| D.34P . Then we use the formula (ID.40p to transform the n-dimensional 
boxes into (n + 2)-dimensional boxes. We find: 



Discs (a'"'''''{- - ++) ) = ieV2) 



,4(12) ]3£ 
[12] (34) 



^Discs (/r'(1324)) + Discs (i^r(1234)) 

^ (7(2,3,4) 



+ E (^Disc3(/2"W) + 4-Discs(Js"(n)) 

,7(1,2)^ ' 



(63) 



3.3.3 Two-cut technique 



We have again two kinds of topologies. The discontinuity, in terms of cut-diagrams, is: 



Disc2 A 



1 scalar / 



+ + 



/ 



pr ^ Di 'pt 




In 



Di Pi 



"1 ,P. 
^1 



+ 



■E 

(7(1,2) 



Pi D'i 

h 

Di 




P{ ^ D \ 
Di + D 
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Here, thanks to the permutations, ah diagrams are doubled. So we gather diagrams and the 
discontinuity becomes: 



Disc2 - ++)^ 

+ (eV2)' V fro (Jl3m+ V 

+ (eV2) 2-. / ^ ^ [i,](^3) + [Ir] {R3) 



[2r](i?4) [2r] I^f (i?4) 



<7(2,4) 



= /i + /2. (65) 

The collection of diagrams, to create gauge invariant trees, has mixed the topologies. We express 
the trees of Ii with the formulas (|F.49| F.50p . For the on-shell trees of I2 , we can use the formula 
()F.56p . but it is not the best way. We obtain directly: 

M(m) ^ [rgil] i {Rq23) i ^ ^ 1 



[r4](i?2) 



^ ' (7(2,4) ^ ^ 3 \ / ^(2,4) 3 

So the discontinuity Disc2 (^A^'^"-'-"'^ { 1"+)) becomes: 

Disc2(^^'^'^''^^(--++)) = {eV2)'Yl E Ti|liy°^''^'^-(^^"(^^^^)) 

(7(1,2)^(3,4) ^ J \ ' 

+ (-^)' E J / ^"Q E ^ ) ([4'^22]^ E ^) '^MM^D 

(7(1,2) V (7(1,3) V V (7(2,4) 3y 

= iVi + iV2. (69) 

Now we simplify the second integral We first introduce two spinors (34) and [21] to build 
the numerator as one product of spinors. We obtain four integrals with the same numerators: 



(68) 



No 



(^^)' T.^fd'-Q ( [3?4l]^ E ^2 ) IM' E ^2 ) ^ m S {Dl) 

(7(1,2) V (7(1,3) V V (7(2,4) 3/ 



(7(1,2) 



(70) 



We have four integrals of rank four. We can use standard reduction techniques to integrate 
them, but it is not the most efficient method. It is better to use the property of the "axis of 
cut" , which is an axis of symmetry. The distribution of helicity is symmetric in relation with 
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this axis. To simplify the expression of the numerator of I2, we use the on-sheh conditions of 
external photons and we note q2 = q- We write the numerator as the product of two equal 
scalars, named P, according to the symmetry of the cut axis: 

(Igi34g321gi 34^321) = {Iq34:q21q34:q21) = (Ig34g21) (1^34^21) = P'^. (71) 

We want to decrease the rank of P and to introduce Df in the numerator. Using gamma matrix 
relations, we obtain: 

(Ig34g21) = 2(g.4)(lg321) - 2(g.3)(lg421) + (/u^ + m^) (13421). (72) 

Now, to continue the simplification of P, we have to know the distribution of photons around the 
loop. The scalar products 2{pj.qi) can be expressed as a sum of denominators and Mandelstam 
variables. Many tensor triangle integrals can be eliminated. A one external mass triangle 
integral with rank one or two can be expressed as a linear combination of one external mass 
scalar triangle and scalar bubbles, according to the formula ()D.32p . However permutations 
allow us to simplify many tensors. The integral (1700 becomes: 



a{l,2) 

/ (I(?34g213g21) - s (^^ + m?) (Ig34g21) -s (^^ + ^n^) {lqUq2l) 
-s {iJ? + m2) {lqUq2l) (Ig34g213g21) -s{^? + m^) {lqUq2l) \ 



(73) 



We apply, again, the development (|72p in each terms of (j73p to reduce the rank of each integral. 
When we have only rank one terms we integrate with the formulas (1D.33| D.34[) . And the 
integral I2 becomes: 

.2,A7(1324).«J,3124„-li^,yM), ,74, 

We gather A^'i and to rebuild the discontinuity Disc2 (j4''™'"''( 1"+)) ^^id we obtain: 

2tu 



Disc2 f^--'-(- - ++)) =(e^/2^^^^^^ ^^^^ 



Disc2 (/4"+2(l324)) + Disc2 (i^4"(1234)) 

<7(2,3,4) 



[12] (34) 

(^DisC2(/2"(n))+4^DisC2(J3'^(7^))) 



a(l,2) 



The reconstruction of the helicity amplitude is easy and we find: 



(75) 



^scalar( , , ^ _^ „-„,2 (12) [34] f 2tu +2/ 



^4 (- - ++) =4 ^« (^-_/^"+^(l324) 
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+ E r-j^i5{u)+^-/uu)) + E ^r(i234) . (76) 

<t(1,2) ^ ^ ,7(2,3,4) / 

This expression is valid to all orders in e. One of the reason of the compactness of the result is 
that we have a symmetry of the helicity structure. Moreover, thanks to the fact that we use a 
four-point scalar integral in n-|-2 dimensions rather a four-point scalar integral in n dimensions, 
we don't have any triangle except the scalar integral J3. 



3.4 Collection of the main result of the four-photon helicity amplitudes in 
massive scalar QED. 

The helicity amplitudes of four-photon scattering are: 
Ar^^^{+ + ++)= Aia' 7t|S^4 (1234) , (77) 



o-(2,3,4) 



y ^ \ I \ ^ ^(2,3,4) 



(78) 



+ E (V^/2"(n) + 4^J^(n)) |. (79) 

a(l,2) ^ 

If we compare the different cut techniques, we see that the four-cut technique is very powerful 
to calculate the coefficients in front the boxes. But we cannot obtain the coefficients in front of 
bubbles and triangles. We can point out that the three-cut technique is sufficient to reconstruct 
all the amplitude. I explain this fact in the subsection 15.31 In the leading order in e, the extra- 
scalar integrals are purely rational (Appendix ID.2p . this is the origin of the rational terms. 
Thanks to the spinor formalism, the results are more compact than all results of four-photon 
amplitude obtained in the past. In the Appendix [Bl I give the massless limit in the leading 
order in e and I find the known results. Therefore, we point out that, in massless theory, only 

the MHV amplitude (^|'^"''^'"( 1"+)) ^ polylogarithm structure in the leading order in 

e. The two helicity amplitudes A^'^"'^"'^ (it -|- +-I-) are only rational terms. 



errrnon 

4 



4 Four-photon helicity ampUtudes in QED: A-^ 

4.1 A{™"(1±, 2+, 3+, 4+) helicity amplitudes 

The two helicity amplitudes A{'^''™°"(1^, 2^, 3^, 4^) in QED, are directly related to the 
scalar QED helicity amplitude A|'^"^"''(l^, 2+, 3"^, 4+) in massless and massive theories: 

Af"™°"(l±,2+,...,A^+) = -2 Ar''^^(l±,2+,...,iV+). (80) 
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This result is true diagram per diagram. To proof this result, we consider a fermion loop with 
N photons entering the loop. We impose, first, that all photons have a positive helicity and the 
same reference vector \R). Therefore we have £ [l---^^], = 0. Now we develop a 

one-loop diagram, called i^/e^"™""^ in QED: 

fermion „A'^ / jn/^ 

(^1 (^1 + "i) •••^TV (^Af + "^)) jv f rnr^^"^ (^i^i-'-^n^n) 

= "W ' = "W ' DI..DI ■ 

(81) 

All terms proportional to m? are proportional to £i.£j = 0, and so vanish. Now if we put the 
explicit formula of the polarisation vectors of each photon ([9]) in (I8ip . then we obtain directly 
that the amplitude of the scalar QED diagram corresponding £)/e'~™o" = _2Dscaiar ^ ^^^1 
diagram and we obtain straightforward ([50]) . Secondly, for — 1_ ++), as we have one 
negative-helicity photon, we impose the reference vector of the positive-helicity photons equal to 
the momentum of the negative-helicity photon. In this case, we have V(?, j) S [1--A^], Si-Sj = 
too, and the proof is the same as the first case. 



4.2 Relation between the QED theories 

We can relate the different QED theories with the Gordon relation. A development of this 
link was initiated in [26]. Currents, in QED, are charged whereas in scalar QED, currents are 
not charged. So to relate the two theories, we have to separate the QED in an uncharged part 
and a charged part, which is the magnetic moment of a gauge field. 

We define the magnetic momenta of a gauge field, with a momentum p and the helicity a, 

by: 

= ea>^^p,e,l = -[f,,^] = -{f^i^-^f^). (82) 

The spinor formulas of the polarisation vectors ([9]), give us: 

= i eV2\p-)ij)+\ and ^/p" = -i eV2\p+){p - \ . (83) 

We consider a vertex in QED between an ingoing photon with a momentum p and two 
fermions with the momenta k and k + p. We decomposed the sum over the two ingoing and 
outgoing currents of fermions in the vertex as the simple vertex of the scalar QED plus another 
vertex called "the magnetic term", just with some gamma matrix relations: 



J[±^, ¥ . , ^ , ^^ ^ 

; ^ff„ — ies„ ^ = —le (2k + p + p) — ^ — ies„ 

—ie 



{k + pf 



{2k + p).e^-^{2^ + ^) + ^^ + /p^) (85) 

((2A: + P) .£p + ^tr^'ep^P.) (86) 

— IP 

(2F + icT'^'^p,)£'^. (87) 



~ {k + pf ' 

In the left hand side of this relation, we have the QED vertex with two currents, which look like 
the currents of fermions. In the right hand side, we recognize the simple scalar QED vertex and 
the magnetic moment of the photon entering the vertex. So we are going to define an effective 
interaction, which describe the QED. 
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We define an effective interaction, described by the vertex Up between a photon, with the 
momentum p and a fermion, with a momentum k: 

Up = -ie{2k'' +pf' + ia'"'pu)ei^ = -ie {2k'' + pf") e^^ + ^p. (88) 



The Gordon relation is written, with this effective interaction: 

[k + pf ' '{k+pf (k + pf ^ 

We are going to show that the relation between QED and scalar QED is complete for a loop of 
fermion with ingoing photons. 

Consider the amplitude 71+72+73+74 ^ in QED theory. If we note Bg^ = 2ie'^rj''^ ei^e2i, 
the double vertex in scalar QED, then the amplitude in QED becomes: 

_ _1 /■ .AtriJWUM , .^j BftrjUM . .^ BftrjUiU^) .^ Bftr {U1U2) 

^ 8 / ^ DlDlDlDl ^ ' DlDlDl ^ ' DIdIdI ^ ' DIdIdI 



a(l.A) 



+ DlDlDl + DlDl + DlDl ' ^'"^ 



To proof this relation , we consider the four-photon QED amplitude: 



fermion _ 

. a n2n2n2n2 • y^^l 



o-(1..4) 12 3 4 



First we begin to develop the amplitude j[f^^"^^°^ [^i 2^ = 16 terms. Then we apply the relation 
linking QED vertex and the effective interaction (j89|) to eliminate all in the numerator 
of each term. In the next step we use again the relation (I89p . which inverts the rotating 
direction in the loop. So, to we find the initial direction, we apply the gamma matrix relation 
: tr (7i....7Ar) = tr (7JV....71). At the end we restore the symmetry to create double vertex B^. 

This result is remarquable. QED is written like scalar QED except for the fact that the 
simple vertex — ie(2/c^+p^) becomes the effective interaction Up = — ie (2A;^ + + icr^^^Pi^). 
This result can be extended to the N-photon one-loop amplitudes. Using this trick, we are going 
to calculate the last helicity amplitude 74{'^^™°"(1~, 2~, 3"*", 4"'"). 



4.3 Af''™''"(l-,2-,3+,4+) helicity amplitude with four-cut technique 

The four-cut technique assumes that all the propagators Dl, Dl, DI^D\ are on-shell. Using 
the formula ([90|) . the discontinuity Disc4 '~~'~'') 

Disc4 - ++)) =-\Y. I d^'Q tr {U,U2UM 5 {dI) 5 {D^) 6 (Z^f) 6 {Dj) . 

(t(2,3,4) 

(92) 



20 



The effective interaction Up = —ie {2k^ + + ia^^pu) £p is a sum of two terms, so the de- 
velopment of the discontinuity gives us 2^ = 16 terms. But, a moment magnetic (I83p is a 
commutator, therefore, a trace of it is zero tr{Mi) = and a trace with two magnetic moments 

with two different helicities is zero too. So the development of Disc4 ^^/<^'"™°"(^ '~~'~^) 

only five terms. The one, with only the QED scalar vertices, is the scalar discontinuity with 
the factor "-2" : 



Disc4 - = -2 Disc4 (^r'"(- - ++ 

+ f d^Q tr (M^M^) ei.qi £2-92 S {Df) 6 {D^^) 6 (Df) 6 {Dj) 

<t(1,2) (t(3,4) 

+ E E -^(-2^e)2i4 f <PQ tr {M^M^) £3.93 £4-94 5 (l??) 5 {dI) 6 {dI) S {Dj) 

a{l,2) a(3,4) 

+ ^ E E -^(-2*e)2i^ /" rq tr (MsM^) e,.qi £2.^3 S {Df) 6 (dI) 5 (Dj) 5 (Z)|) 

ct(1,2)ct(3,4) 

+ IY. Yl -^(-2^6)2^4 j cTQ tr iM,M2) £3.92 £4.94 5 {Df) 6 {dI) 5 {Dfj 5 (Z)|) 



a{l,2) <7{3,4) 

1 scalar / 



= - 2 Disc4 - +h + l2 + h + h. (93) 

The factor "2" in front of the scalar discontinuity, comes from the fact that we need two complex 
scalars to build a fermion and the sign "-" comes from the fact that we change a fermion-loop 
into a boson-loop. We find this factor in the supersymmetric decomposition ([T]). Now we have 
just to calculate all trees containing in ([93]) to obtain the discontinuity. We, first, compute the 
traces of magnetic moments using the definition ()83p and then we simplify the formula of trees 
in dMD with (IR45]) . We obtain: 

Disc4 - ++)) = - 2 Disc4 - ++ 

■ (12) [34] 
[12] (34) 



+ (^^^'ttIM^ I E (Jr(1234)) - 2Disc4 (ir+'(1324)) 1 . 



^(2,3,4) 

(94) 



Now we use the two-cut technique to calculate the same amplitude. 

4.4 2 ,3+, 4+) helicity amplitude with the two-cut technique 

The QED discontinuity Disc2 '~~'~^) 

C/2f/3f/4)w^2N W^2^ , ■4tr{UiU2UsU4) ,.^2^ .rr.2^ 

Bisc2\^Ai (--++)j=--2^ dQi -2-2 6{D2)6{D^)+i 6 [D-^) 6 [Dr^) 

(7(2,3,4)-^ 13 2 4 



D2 1 



(95) 
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The computation give us directly: 

Disc2 _ ++)) = _ 2 Discs [Ar'^'i- - ++; 



+ I E (Jr(1234)) - 2Disc2 (ir^(1324)) 



^o-(2.3,4) 

(96) 



4.5 Collection of the main result of the four-photon helicity amplitudes in 
massive QED. 

With the two-cut technique or the four-cut technique, we find the same discontinuity of the 
MHV four-photon amphtude in massive QED. The reconstruction gives us: 



++) = -2^r'"'(— ++)+4ia'|^|^s I J2 '/r(1234)- 2 /4"+2(i324)j . 

(97) 



^(t(2,3,4) 

If we take the formula of scalar amplitude already calculated (|75|) . the QED amplitude is : 



Af~i--++) = -8^«^j||||(|^/r^(1324)+ Y: (V^/2"(n)+4^J3"(n) 



+ [K2{1234) - - Jr(1234)) . (98) 

a(2,3,4) / 

The known massless limit is given in the Appendix [Bj 

5 Discussions 

5.1 On the analytical structures 

At one loop order, the helicity amplitudes of the four-photon process have the same structure 
for QED or scalar QED theories. The first two are only a rational term whereas the MHV 
amplitude has a polylogarithm structure, carried by the scalar integrals I^'^'^ and I^- 

Now we are going to prove that in massless QED and in massless scalar QED, and therefore 
in massless supersymmetric QED-'^^^, four-photon amplitudes could be decomposed without 
triangle. It comes from the fact in the decomposition ([2]), the infrared divergence are carried 
out by scalar triangles. Consider one diagram, named L of a four-photon amplitude in QED (we 
have exactly same proof for the scalar QED). The numerator of fermion's propagators implies 
that all IR singularities vanish. Now consider a sub-diagram of L by removing propagators. 
We have four sub-diagrams. After removing a propagator that connects two photons with 
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on-shell momenta pi and p2, the amplitude depends only on the sum of the two momenta 
g = Pi + P2 with non zero ^ 0. The two photons around the pinched propagators have the 
behavior of one massive photon. However, the mass of a massive entering particle regularizes 
IR divergences. Therefore all sub-diagrams of L are not IR divergent. So as each sub-diagram 
is a tensor triangle and doesn't have any IR divergences therefore each sub-diagram cannot be 
decomposed with scalar triangle. Finally we can not have any triangle in massless theories. But 
in massive theory, this argument is not exact, because, the infrared divergences exist only in 
massless theories. However no triangle is expected, except some extra-scalar triangles with the 
extra dimension term fj? + m^. 

The bubbles have two origins. The first origin is the decomposition of three-point tensors 
integrals and the second origin is the UV divergences of the loop. Here we have no triangle, so 
the contribution of bubbles coming from the reduction of triangles is zero. Now we study the 
UV limit of one diagram in scalar QED. We use the notations introduce in subsection 12. 3[ For 
example, we have: 



The tensor integral is UV divergent, therefore the reduction creates bubbles. But, whatever 
the helicity amplitudes, the contraction of tensors ei £2^364 [ij^'^ij^'^ + ■q^Pr]'^'^ -|- rqt^^rfP^ is zero. 
So we conclude that each diagram of the four-photon amplitude is UV finite, thanks to the 
gauge invariant. The UV finiteness express by compensations of the divergences of the bubbles. 
We clearly observe this phenomenon in the MHV amplitudes. But, for the first two helicity 
amplitude ^ (it + we have at least three positive-helicity photon, so all the discontinuities, 
in four dimensions, are zero and it implies that the coefficients in front of each bubble for each 
diagram is zero. 

5.2 On the origin of the rational terms 

In the very simple example of the four-photon amplitudes, we point out that the rational 
terms come from the extra-dimension integrals Jj and Ki. We can discuss the origin of those 
integrals. Consider a one loop diagram of the four photon amplitude in the four-dimensional 
helicity scheme, described in the paragraph 12. 3[ We can write a diagram as: 



where Num(Q) = ei.Qie2-Q2£'i-Q?,£i-QA- However the regularization scheme imposes that the 
vertices are in four dimensions. And the 4-dimensional space and the — 2e-dimensional space 
are orthogonal therefore, £i.Qi = £i-qi- The numerator is actually a function of the four 
dimensional part of the loop momentum: Num((5) Num(g). Moreover the denominator of a 




(100) 



(99) 



(101) 




(102) 
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propagator is spelt: Df = qf — — m?. During the reduction some squares of momenta appear, 
like g?, in the numerator. To rebuild a denominator, we subtract and add the mass: m? + ^^: 

The integrals with /i-^ are only rational, this is the origin of the rational terms. In this case, it 
is very simple to find the rational terms. We have just to shift the mass of the scalar: 

m^ + ^?. (104) 

The four-photon amplitudes are a special case where the particle in the loop is a scalar or 
a fermion. Consider the case where we have a photon propagator in the loop. The internal 
photon is in n dimensions and its propagator, proportional to the metric r}^'^ . The contraction 
of this metric in n dimensions with loop momenta Q/^ creates some //^ terms because there are 
not enough vertices to reduce all the propagators in four dimensions. In this case, we cannot 
associate the /U^ terms with a mass. 

The four-photon amplitudes have no ultraviolet and infrared divergences so the amplitudes 
don't depend on the scheme of regularization. In the case where the amplitude have diver- 
gences, which are not regularized, then the rational terms depend on the scheme. In [25], we 
have methods to change the scheme of regularization without much eff^ort, just by adding or 
subtracting a rational term. 

5.3 On the multi-cut techniques 

In this work we apply three kind of unitarity-cut techniques with two, three or four cut 
propagators. We first note that, the more there are cuts, the more we have on-shell conditions 
and the simpler is the computation. But the more we have cuts, less coefficients in front of 
the scalar integrals could be calculated. Actually the four-cut technique is very powerful to 
calculate the coefficient in front of the four-point scalar integrals. The three-cut technique is 
sufficient to calculate the coefficient in front of scalar triangles, scalar boxes and scalar bubbles. 

However, the fact that we can calculate the coefficient in front of bubbles with three-cut 
technique, is a peculiarity of the four-on-shell-photon amplitudes. Consider one four-photon 
diagram which is made by four photons ingoing into a fermion loop. We assume that we apply 
the two-cut technique to it. We call the axis of cut, the main cut. This axis shares the four 
photons in two trees of two photons. Each group is the same invariant. Then we add a cut. So 
we cut one tree with two photons into two trees with one photon. As one on-shell photon cannot 
constitute an invariant, therefore when we cut, we don't divide an invariant into two invariants. 
So it remains only one invariant, which is the one of the two cut technique. Therefore, we don't 
touch the analytic information contents in the branch cut and don't loose information when we 
extend the two-cut technique to the three cut technique. This fact explains why the two-point 
functions, which respect the main cut, are not spurious in the three-cut technique. 

6 Supersymmetric amplitude A^^^ 
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We use the super symmetric decomposition ([T]) to extract directly the supersymmetric am- 
phtude A^=^. Since the j\f^^^^°^ obeys to the supersymmetric decomposition ([T|), that means 
that we can identify the A^^^ , without computing all diagrams. So, with the formula ()80p . we 



can identify the ^(i + ++) and, with the formula (j97p . we can identify A^=^{- - ++). 
We obtain: 

Af=i(+ + ++) = 0, (105) 

A^=^{- + ++) = 0, (106) 

A^=\- - ++) = 4 i «'|^|^^ 1^ E mm - 2/r'(i324) j . (lo?) 

In massless case, in the limit e — > 0, we have: 

A^=\+ + ++) = 0, (108) 

A^=^{- + ++) = 0, (109) 

A^=^{--++) = -8 i o? rloilofi g^r^^(1324) + 0(6). (110) 



There is full agreement with [5]. With a massless or massive loop momentum, the supersym- 
metric amplitudes have no rational term, no bubble, and no triangle, only boxes. 

We are going to prove that diagrams of the four-photon amplitudes in supersymmetric QED 
: M = \ are UV finite. We identify the decomposition of a fermion loop (I90p and the formula of 
the supersymmetric decomposition ([TJ. As the trace of one magnetic moment is zero, therefore 
we see that all the terms belonging to the supersymmetric amplitude have at least two magnetic 
moments. So we can do the power counting of one of those terms. We define r the power of the 
loop momentum of the A^-photons amplitude. We have: 

r = n-\ + N = l-2e-7V. (Ill) 

So, if > 3, therefore the loop is not UV divergent and so diagrams of the four-photon 
amplitude in supersymmetric QED-'^^^ are UV finite. 

In the last section, we show that the diagrams in QED^'^^^ have no IR divergence and 
therefore no triangle. The bubbles in the supersymmetric amplitude could come only from the 
UV structure. But we see that each diagram has no UV divergence, so they are no bubble in 
each diagram. We can observe it with standard reduction of the four-photon amplitude. There 
are some interferences in the loop between bosons and fermions, which reduce the UV power and 
eliminate all bubbles. The interferences create magnetic moments, which are gauge invariants. 
Interferences increase the power of the gauge invariance. 

In the next section, we calculate, the most simple helicity configuration of six-photon am- 
plitude in massive theories. 



7 The first helicity ampUtude ^6(+ + + + ++) 
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In |27l EH] ! the six-photon helicity amplitudes in massless and massive loop were numerically 
computed. Here we obtain an analytic expression of the most simple helicity amplitude, all the 
six photons have a positive helicity for a massive loop. 



A six-photon one- loop diagram is not IR/UV divergent, so in this part, the dimensional 
regularization is not needed here and the integrals are in four dimensions. With standard 
techniques, we show that this amplitude has neither bubble, nor rational term and nor triangle 
with one and two external mass. We verify it explicitly by the computation. 



Thanks to the supersymmetric decomposition ()80p . the scalar amplitude gives us directly the 
fermionic amplitude and the supersymmetric amplitude. Now consider a diagram and we apply 
the two-cut method. There are two kinds of discontinuity. The first kind of discontinuities sep- 
arate the six photons in two groups of three photons, whereas the second kind of discontinuities 
separate the six photons in a group with four photons and a group with two photons. We don't 
have bubble contributions so the second kind of discontinuity is better because with only one 
cut we can have the coefficient in front of all kind of the scalar integrals. The problem of the 
first kind of discontinuities is that we cannot have the coefficient in front of triangle with three 
external mass. Let us cut the diagrams in the channel s^q: 

Disc2,s56 (^r^'^'' (+ + + + ++)) = i-2iefi^ J A^tree(l+,2+,3+,4+)^t„e (5+ , 6+) <5 I?!) 

(112) 

In this part, to reduce, we note 6 (Z)|,Z)|) = 6 6 We need on-shell trees with two 

photons and four photons. The computation give us: 

A,ree (5+,6+) = V fll^^ = -m' ^ V i., (113) 
^..4l-,2-,3-,4-)= '"^^1?^^^?'^"'^^ (114) 

cr(l, 2,3,4) 1 2 3 

[12][34] 1 , 4 [4(4 + 3 + 2 + l)gol] 1 



4 

m 



V ^ ^ - +m^ V 

^ 4(12)(34) D2^2^2 + ^ -1-2-3 

(7(1,2,3,4) \ ' 12 3 a(l, 2,3,4) \ /\ /\ / 12 3 

(115) 



We can find some quivalent formulae in [29]. We put those trees (I113| 115]) in the amplitude 
(|112|) . We have to integrate tensor hexagons rank one, which give only four point scalar integrals. 
So we are going to calculate the coefficient in front of each scalar box. We begin to calculate 
the discontinuity of the two adjacent mass four point scalar function : 

Cham 2 Q3 ^« 

P2 \ ^ / P5 




where we have two invariants si2, sse : 

Disc4,si2,.56 (^r'"') = j d\ Atree (l+,2+) \{e,.qi Atree (5+,6+) 5 {dI Dj, Dj, dI) 



i=3 



(116) 
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[12] [34] [56] 
(12) (34) (56) 



Disc, 



4,S12,S56 



(117) 



We use the trees (|113p . the computation is straightforward and the reconstruction too. Ah 
information of the two adjacent mass four point functions is hold in the hexagon. Now, we 
compute the coefficient in front of the one mass four point function: 




chain 1 



The discontinuity is: 



Disc 



scalar 



J A Aarbre (l + , 2+ , 3+) £+ .^44 -954 -96 ' ^1' ^1) • (US) 



To not break the symmetry, we use the four-photon scalar trees (jll5p . The trees with a even 
number of photons are simpler than the trees with a odd number of photons. We simplify 
thanks to the permutations and after reductions, we obtain: 

Disc4,.,,3 (^6) = - {eV2fm^ n^lf^Tl E Disc4,.,,3 i^e 

(t(1, 2,3,4) 



(e\/2)"m 



(12)(34)(56) ^ 

^ ' ct(1, 2,3,4) 

[56] V- 



4 L^"J \ " 

(56) ^ 

^ ' <t{1,2,3,4) 



(12)(23)(34) Df 



HDlDlDlDl), 
(119) 



where M = pi + p2 + P3 is the momentum of the external mass. Rotating the gammas matrix, 
and using the on-shell conditions, the numerator becomes: 



(54^56165) = L>f (54^65) +m2(54165) - (54^^3165). 
So the discontinuity is : 



(120) 



Disc 



4,si23 1^6 



A scalar \ 



J 



(ex/2)* 



m 



[12] [34] [56] 
(12)(34)(56) 



Disc, 



4,si23 (-^6 



+ {eV2ym 



(54^/65) 



-DisC4,si23 {III (^123)) • (121) 



2(12) (23) (34) (45) (56) (61) 
Finally the compute the discontinuity of the two opposite mass four point scalar integral: 

P6 

Pi 
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and the discontinuity is: 

DisC4,si2,S45 (^r'"") = / d\ Atree (l + , 2+) £+ Qs^tree (4+, 5+) S+gg S DI DI D^) . 

(122) 

To not break the symmetry of the scalar integral we use the three photon scalar trees: 
We apply twice this formula and, thanks to the permutations, the discontinuity gives: 

Disc4,.,„.« (^r'"'^) 

cr(l,2,3)cr(4,5,6) 

V2fni'' y [ d\ , J.^^f^^^S^/^^^^, , 5 iDlDlDlDl) . (124) 
^ ^ J ^(12)(23)(45)(56)L>?i:>| ^ ^ ^ ^' ^ ^ 



o-(1,2,3)<t{4,5,6) 

We simplify the numerator by rotating the gammas matrix: 



f3^211[604l - [3216543] , , [3(1 + 2)5(5 + 4)3] 
2 [36543] ^2 [3(1 + 2)^26543] 

+ ^1'^^ MW " ' — WW) — ' ^ ^ 

and we obtain directly: 

T~»;r,^ / Ascalar \ - 1^ /^^6_6 [12] [34] [56] 

JJisc4,s-^2,s45 j — - yey A) m ^-|^2)(34)(56) ^•**i2'''45 V-'e ) 

, / A^n6 4 [3(1 + 2)6(5 + 4)3] ^. , 
+ iey2Yirr , \ ; — , , ' — , , ' , — -D1SC4 s,, (/roR(si2, S45)) . 
^ ^ 2(12) (23) (34) (45) (56) (61) 4,^12,54. v 4,2BV i^. 457; 

(126) 

In each discontinuity, we have the trace of the scalar hexagon. The entire reconstruction is 
immediate: 

.scalar ^ _ {eV^frif_ [12] [34] [56] 

96.^ .a,2tl,.a)^^2)^^^)^^6) ' 
W2)6m4 [3(1 + 2)6(5 + 4)3] ^„ 

96vr2 ^^^^Z.^^^^^^2(12)(23)(34)(45)(56)(61)'^^^^ 

W2)^ ,p (54^65) 

96vr2 2(12)(23)(34)(45)(56)(61) ^ ^ 

a-(l,2,3,4,5,6) 

The coefficient in front of each four point scalar integral could be written as: 

det (S) 

(12)(23)(34)(45)(56)(61) ' 

where S is the kinematical matrix of the massless scalar integrals. If we put all the fonctions 
in n + 2 dimensions, therefore all three points functions vanish. Thanks to the supersymmetric 
decomposition ([50]) . we have: 

+ + + ++) = -2^^^'^''^'^ (+ + + + ++), (129) 
A-^=\+ + + + ++). = (130) 



(128) 
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8 Conclusion 



In this paper, we have calculated all the four-photon helicity amplitudes in massive and 
massless QED, scalar QED and supersymmetric QED^'^^^. To compute them, we use two 
very powerful methods: the unitarity-cuts and helicity amplitudes accompanied with the spinor 
formalism. So we don't need any standard reduction method. Thanks to four cuts in four 
dimensions, we obtain the coefficients in front of boxes, thanks to three cuts we obtain the 
coefficients in front of triangles and bubbles. The extension in 4 — 2e dimensions of the unitarity- 
cuts allows us to calculate straightforward the rational terms. 

In this example, we have simplifications because we have only four massless external legs. 
So we don't need the two cuts methods to compute the coefficients in front of bubbles. 

Moreover, the transformation of the four point scalar integrals in n dimensions in four point 
scalar integrals in n-|-2 dimensions reduces strongly the final result. Indeed this transformation 
allows us to separate the infrared structure of the amplitude. This separation generates many 
explicite compensations. 

As we have some very compact expressions of the six-photon amplitudes in the massless 
theories [33], we hope to obtain expressions in massive theories, by understanding of the origin 
of the rational terms. Thanks to the two-cut techniques, we could calculate the first of the 
four six-photon helicity amplitudes. In a next paper, we develop the calculation of the next 
six-photon helicity amplitudes. 
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We give, for sake of completeness, the vertices in QED and scalar QED, and then the known 
results on the four-photon amplitudes are given in the second appendix. The third appendix 
recall the reduction of tensor integrals and extra-dimension scalar integrals. Moreover, in the 
fourth appendix, we give the definition of the master integrals used in this paper is recalled and 
in the next appendix the reduction of the pentagon. Finally, we give the proof of the amplitude 
of each chain of photons used in this paper. 



A Vertices 

The QED vertex is: 



APPENDIX 




(A.l) 
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whereas the two scalar QED vertices are: 

' \^ ^ =-ie{k + {p + k)Y ^-ff" =2ie2ry^^ (A.2) 

k + p/ ^ 



B Massless limit of the four photons ampUtudes 



We take the massless limit in the leading order in e of the results in the section 13.41 We use 



the Schouten Identity, and we find that (34) ' (34) (231) invariant by permutation. So the 
helicity amplitudes are: 



Ar'-(- - ++) jilliy {1 - '^-fn^^^m + ^ m^) - m))] + o{e\ 

(B.5) 

4™"(± + ++) > -2Ar^'"'{± + ++), (B.6) 

- ++) jilliy {1 + ^^^r^'(i324) + ^ (/^"(n) - im)] +o{e). 

(B.7) 

We find the results known in the massless limit [U \E\ . 



C Definition of the master integrals 



In this appendix, we give the definition of the master integrals used in this paper. We write 
Gi the Gram determinant and Si the kinematical S-matrix. Consider a scalar integral: 

pi 



■71 

P3 



We define the Gram and kinematical S-matrix by: 

G^j = 2pi.p,, (C.8) 
Sij = {qj - Qif - mf - ruj , (C.9) 

and we note the spatial integral: 

Finally, in all formulae, the analytic prescriptions are: 

s — > s + rr? — iX. (Cll) 
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C.l Two-point functions 



In a massless theory, the two-point scalar integral in n dimensions is: 

l2is) = J^^^^i-s)-' = ^-ln(-s) + 2 + 0(e). (C.12) 

and in n + 2 dimensions, the two-point scalar integral is: 

/"+2(s) = !:L -(-s)^-\ (C.13) 

2 ^ ^ 2e(l-2e)(3-2e) ^ ^ ^ ' 

In a massive theory, the two-point scalar integral in n dimensions, in the leading order in e is: 

I^is) = + 2 + / (^) + 0(e). (C.14) 

where p = 1 — . Thanks to the small imaginary part: s — > s + iA, — > m? — iX and > 0, 
we have: 

if s < then / (^) = ^ln( ^~| )-iA (C.15) 



if s > W then / (^) = ^ (^In (^^^) " (C-16) 
if s G [0, 4m2] then / (^) = -i^/^ In ( ] ] . (C.17) 



The determinants are given by: 



det(S2) = s (s-W) , (C.18) 
det(G2) = s. (C.19) 



Most of those results comes from ^ |6j . 

C.2 One external mass three-point functions 



Pi 

In a massless theory, the one external mass scalar triangle in n dimensions is: 

4" (») = 5^ = i - 1" (-) + ^) + o (.) , (C.20) 

and in n -|- 2 dimensions, this triangle is: 
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In a massive theory, the one external mass scalar triangle in n-dimensions, in the leading order 
in e is: 



where p = \/ 1 — . The determinants are given by: 



det(53) = 2s'^m^, (C.23) 
det(G3) = -s^. (C.24) 



C.3 No external mass scalar four-point function 



pi 

t 



S = Si2 
t = Si4 
U = Si3 



P3 

In a massless theory, the no-external-mass scalar box in n dimensions is: 



^4(«>*) = ^5{(-^)"^ + (-*)"1-^^o(5,i), (C.25) 



where F(){s,t) = — |ln^ + In n + 2 dimensions, this box is: 



In a massive theory, the no-external-mass scalar box in n dimensions, in the leading order in e, 
is: 

where: 

H{X, Y) = (in (l - ^) In ( - — 



y-x+J \x--yj \x+-yj \y-x 



Li2 ( ^— ]-Li2( ) + Li2 ( ) \ , (C.28) 



with x± = — ^1 lb \/l — 4Xj and ^ ~ 2 ~'~ ~ ' '^'^^ standard definition of the dilog- 
arithm is Li2{x) = — (it in(i^ xt) ^ rpj^^ determinants are given by: 

det(54) = ts[-Arr?{t + s) + ts\, (C.29) 
det (G4) = -2st(s + t) = 2stu. (C.30) 



D Reduction of integrals 
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D.l Reduction of tensors integrals 



In this appendix, we give the reduction of the hnear-tensor two-point integrals, hnear-tensor 
one-mass three-point integrals and linear-tensor no-mass four-point integrals. In the argument 
of each integral, we give the numerator of the tensor integrals. We can find techniques of 
reduction in [31] . 

■ D CX"' >C< 

Pi ^ ?1 >2 Pi Pi 

Figure 2: Kinematics of bubbles, one external mass triangle and no external mass box. 



(D.31) 
(D.32) 



I2 {4) = 






^« - + 




-l-itl2{1234) 




-i^(tJr(1234) 



^^^^/r(1234), 

(D.33) 

- i^l±^ J«(1234). 

(D.34) 



D.2 Reduction of extra-dimension scalar integrals 

Here we give several formulas to calculate extra-dimension scalar integrals. Most of those 
results come from [3] . We can relate the extra scalar integral and the scalar integral in n + 2t 
dimensions: 

Using this formula we can calculate easily Jm and function. Here we give some formula 
useful for this paper: 

(_e)/^"+2 = + O(e), (D.36) 
(-e)(l-6)/r' = -\ + 0{e). (D.37) 

We have: 

= m'l'}, + {-e)P/\ (D.38) 
= m^r^-2n?tr+^ + {-e){l-e)Il+\ (D.39) 

And the following relation holds between the J4 and I^^"^ functions: 

J4"(1234) = -£/r(1234) + ^I2{s) + ^lUt) - ^/r'(1234). (D.40) 
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E Reduction of the pentagons and scalar hexagons 



The method of reduction comes from [31]. The expHcit reduction of massless scalar pentagon 
and massless scalar hexagon is in [TJ [30l [32] . 

We note /g the hexagon in n dimensions and (sij) the n-dimensional one-mass scalar pen- 
tagon. This pentagon is obtain by removing the propagator, of the scalar hexagon, between the 
external momentum pi and pj. We note Sqij the kinematical matrix of the hexagon and S^ij 
the kinematical matrix of the pentagon. 



According to [30l[3T], the exact reduction of the hexagon is: 

6 



-^6 = X] ^^kl^^5 (si l+l) , 



k,l=l 



where the kinematical matrix is: 



56 








S23 


S234 


S16 


\ 




A 


1 


1 


1 


1 


i\ 











S34 


S345 


S12 




1 


1 


1 


1 


1 


1 


S23 











S45 


■5456 




1 


1 


1 


1 


1 


1 


S23A 


S34 











.556 


1 


1 


1 


1 


1 


1 




■S345 


S45 













1 


1 


1 


1 


1 


1 


\ 


S12 


■S456 


S56 





) 




U 


1 


1 


1 


1 


1/ 



kl 



The reduction, in the leading order in e, of the pentagon is: 

5 



kl=l 



(E.41) 



(E.42) 



(E.43) 



Where I4 (s, si is the four-point scalar integral, obtain by removing the propagateur between 
the legs with the momentum pi and of the pentagon I^{s). The kinematical matrix for the 
pentagon 15(512) is: 



/ S34 S345 Si2\ 

S45 S456 

S34 S56 

S345 S45 

\S12 S456 S56 / 



2m^ 



/I 1 1 1 1\ 

11111 
11111 
11111 
Vi 1 1 1 1/ 



(E.44) 



For another pentagon, we permute the labels. The hexagons and pentagons have no infrared 
divergences. We keep only the "finite parts" of box functions. This part is usually called 
where the subscripts i represents the number of external legs. We found those definitions in [7]. 



F Computation of the on-shell trees using in the paper 
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F.l On-shell trees with two positive- helicity or two negative- helicity photons 



Proposition 1: 

Consider two positive- helicity photons with the momentum pi and p2, join and connect them 
by an on-shell propagator, therefore we have: 

(K5j2) _ [12] 

{Rl) {R2) ' (12) ^ ' 

In the case of two negative-helicity photons, we have: 

[IR] [2R] ' [12] ^ ' 

Proof : We assume to have two photons with a positive helicity. We note qo = qi — pi- And 
all propagators are on-shell, therefore 2(pi.qo) = 2(p2-92) = 0. Using this trick the amplitude is 
spelt: 

{Rqil) {Rq22) ^ {Rqol2q2R) ^ 2{p2.q2){RqolR) - 2ipi.qo){Rq22R) + {Rlqm'^R) 
{Rl) {R2) {RI){R2){12) {RI){R2){12) 

(F.47) 

qfjRUR) {f.^ + m^){R12R) [12] 
{R1){R2){12) {R1){R2){12) '{12)' 

There is the same demonstration with negative- helicities photons. 
Proposition 2: 

Consider a chain of two positive-helicity photons with the momentum pi and p2 surround with 
two on-shell propagators, therefore, we have: 

If we have two negative-helicity photons, the chain is: 

If the joined propagator is put on-shell therefore we find formula ()F.45| F.46p . 
Proof : 

^ {Rqil) i {Rq22) ^ ^ i {Rqol2q2R) 

.^2) ^^^^ ^^2) ^^^^ {Rl){R2) {12) Dl ^ ■ ^ 

i 2{p2.q2){RqalR) - 2{pi.qo){Rq22R) + {Rlqoq22R) 



E 



(m)(i?2)(12) Dl 



E 



(F.52) 

- Dl) {RqolR) - {Dl - Dl) {Rq22R) + ql{R12R) 



^^^^^^ {R1){R2){12) Dl 



(F.53) 
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^^^^ (m)(i?2)(12) Df ^ ■ ' 

^ ^ -(1,2) ^1 

Most of the reduction come from the permutation. For photons with a negative hehcity, the 
proof is the same. 



F.2 On-shell trees with one positive-helicity photon and one negative-hehcity 
photon 



Proposition 3: 

Consider a diagram with a scalar hne. On this scalar line with have first a chain with one 
on-shell negative-helicity photon and one on-shell positive-helicity photon. We assume that the 
momentum of negative-helicity photon (respectively positive-helicity photon) is: pi (respectively 
P2)- Moreover we assume that a third on-shell photon, with the momentum p^, is ingoing in the 
scalar line just after this chain (fig. [3]). The propagators around the chain are on-shell Dq = 
and D2 = 0. The amplitude of this diagram is: 



Pi P2 P3 




Figure 3: Chain composed with two photons, with two different helicity, following with one photon. 



r2]{Rl) ^ [rqil] i {Rq22) _ i / (1^2232) - I)| (1^22) + (/i^ + m^) [231] \ 



2(i?2)[lr]^^^^^ [Ir] Df {R2) (231) \ Df J 

i_ f -Dl{lq22) + (lg2312) + {^P+m'^) [231] \ 

(231) )' 

(F.56) 

Proof : We choose \R) = |1) and |r) = |2). The computation gives us: 

[rgil] i {Rq22) ^ _z(l^ V — = ' ( (^^^231^22) (Ig2231g22) ^ 



(F.57) 



(F.58) 



2(pi-g2)(lg2232) - 2(p3.g2)(lg2212) + 2(p2.g2)(lg2312) - gi(12312) 
512(231) V Dl 

i ( 2(R.g2)(lg2232) - 2(p3-g2)(lg2212) + 2(p2-g2)(lg2312) - gi(12312) 
512(231) 1^ / 

i f (lg2232) - 1^1(1922) + (/i^ + m^) [231] \ i_ f -D|(lg22) + (1^2312) + (^^ + m^) [231] ' 

(231) Dl )~ (231) 

(F.59) 
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G Four-cut technique for ^|ca/ar^ — ^ 



In this appendix, we want to prove that, the integral (j42p could be calculate without inte- 
gration formulas, like (lD.33p . We use just the four on-shell conditions. We want to calculate 
the tensor integral: 

(^^)' E / d-Q{lq2232) + m') 5 {dI) 5 [dI) 6 (Z^f) 5 {Dj) (G.60) 



<7(2,3,4) 



corresponding to the graph: Di - r j -Dl 



The four-cut technique [151 lEl ttZ] says that the integrals (IG.60p is 

,2 I ^2 



where are the solutions of: 

-D? = 4^ (52 - P2)^ - {m^ + Ai^) = (G.63) 

1)2 = ql-{m^ + /i^) = (G.64) 

Dl = Q ^ {q2+ pzf - (m^ + /i^) = (G.65) 

-D| = ^ {q2+P3+ PA? - [m^ + Ai^) = 0. (G.66) 



To solve this system of equations, we choose a basis of the four-dimension Minkowski space: 
B = {p2,p^,{2^'^3),{3^'^2)}. In our case, qf^- is a four-dimension vector, therefore, we can 
project it on the base B: 

= a,pi^ + h P^s + I (27^3) + |(37^2). (G.67) 

So to know the vector qt^-, we have to calculate, the four coefficients ai,bi,Ci and di. The 
conditions ([06311 and (fClM]) impose: 

(Q2i - P2? = m^ + f,^^2 {p2.q2i) = 4^ 6i = 0. (G.68) 
The conditions (fClM]) and (IGTHKI) impose: 

(g2i + P^f = m^ + fi^^2 {p2.q2^) = ^ a, = 0. (G.69) 
The second condition ()G.64p imposes: 

q2i = + /i^ 4^ ad, = —. (G.70) 
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And finally the two conditions (|G.65p and (|G.66p impose: 

{q2i +P3+ Pif = m^ + fi^^ q(243) + d,(342) = -S34. (G.71) 

So we have = 6j = and the two equations ()G.70| G.71|) gives us q and dj. q2i is totally 
define. Now we insert the decomposition ()G.67p in the equation (|G.62p and we have: 

(1 ^ 92.232) = -t(123) ^ c. = st^^. (G.72) 

i i 

We input the last result in the integral ()G.62p . and we obtain: 

(^^)' E I ^'^^^1^2232) {y? + m') 6 {dI) S {dI) 5 (Df) 6 {Df) 

ct(2,3,4) ^ ' 
<t(2,3,4) ^ ' 



This result is the result which we have obtain with the classical integration (j43p . 
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